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ABSTRACT 



We develop a new and general method to calculate the effects of CP violation 
from extensions of the standard model on the mechanism of electroweak baryoge- 
nesis. We illustrate its applicability in the framework of two-higgs doublet models. 



1. Introduction 



It has recently been convincingly established that electroweak baryogenesis 
due to the mixing of three generations of quarks in the minimal standard model 

is unable to account for today's baryon asymmetry!^'^'''' This result is the direct 
consequence of a sharp conflict between the rapid quark-gluon interactions and 
the far too slow processes of quantum interference through which the phase of the 
Kobayashi-Maskawa matrix can emerge into the physical world!*' 

This result has long been anticipated and a plethora of new sources of CP 
violation have already been contemplated. There is a clear contrast between these 
new sources of CP violation and the standard model CKM phase. The latter is only 
physical because of the charge current interactions. In contrast, in many models 
with new sources of CP violation, during the weak phase transition some mass 
matrix has a space dependent phase which cannot be removed since making the 
masses real and diagonal at two adjacent points x and x-\- dx requires, in a space- 
varying background, two different unitary rotations Ux and Ux+dx- The relative 
rotation U~^Ux+dx yields a new interaction which can generate a CP violating 
observable. Because large physical CP violating interference effects can appear in 
the phase boundary where the particle masses are space dependent, they play a 
dominant role in the mechanism of electroweak baryogenesis. 

At the electroweak phase transition, non equilibrium CP violating effects are 
largest inside the wall of a bubble of broken phase expanding in a thermally equi- 
librated plasma. Two apparently distinct mechanisms of electroweak baryogenesis 
have been proposed in the hteraturej^'the so-called "thin wall" and "thick wall" 
scenarios. They are characterized by the conditions £/L^ 1 and i/L respec- 
tively, where L and i are the wall thickness and a typical mean free path of the 
particles relevant to the scenario. 

In the "thin wall" scenarios, coefficients of reflection and transmission are com- 
puted. For those species whose interactions with the wall take place through a 
complex mass matrix, these coefficients assume different values for a particle and 
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its CP conjugate. The resulting asymmetry is then convoluted with an incoming 
thermal flux to yield a CP violating source which moves along with the interface. 
This method has been applied to exotic fermions with a majorana mass and to the 
top quark in two-higgs doublet models!*^' 

In the "thick wall" scenarios, thermally- averaged local operators have been 
written which couple the baryon current, or a related current, to the space-time 
derivative of the mass terms. These operators act as CP violating sources deflned 
at every point of the interface. These are the scenarios proposed in the context of 
the two higgs'^'* ''' and supersymmetric models!"'"' 

It has recently been understood that in both scenarios these sources are to 
be inserted in a set of coupled rate equations which allow for the CP violating 
charges to be transported elsewhere in the plasma. Transport greatly enhances the 
final baryon asymmetry since anomalous electroweak baryon violating processes 
are suppressed in the wall and in the broken phase but are relatively rapid in the 
symmetric phase!^^'^^'"' 

This dichotomy in the formulation of electroweak baryogenesis mechanisms 
obviously reflects limitations of the computational techniques and does not do 
justice to the underlying physics of CP violation. It is unsatisfactory for a few 
reasons. 

(a) In practice, the thickness of the wall is neither very small nor very large com- 
pared to a mean free path. The "thin wall" approximation can overestimate 
the magnitude of the actual baryon asymmetry produced while the "thick 
wall" approximation can underestimate it. 

(6) An interpolation between these two limits is required as the CP violating 
sources are the only inhomogeneous terms in the rate equations so that the 
uncertainty in the baryon asymmetry is directly proportional to the uncer- 
tainty in their determination. 
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(c) For any given energy, particles moving at an oblique enough angle relative to 
the wall are likely to be scattered while inside the wall. Hence the integral 
over all particle momenta will include both particles which scatter inside the 
wall many times and those which do not scatter at all and the "thin wall" 
limit is never fully applicable. 

(o?) It has previously been argued'^ that for a sufficiently thick wall it is a 
good approximation to compute the effects of the nonequilibrium CP vio- 
lation by adding time varying CP and CPT violating terms to an effective 
Hamiltonian which is treated as approximately spatially constant, and by as- 
suming that the local particle distributions relax towards thermal equilibrium 
with this Hamiltonian according to some classical rate equations. However, 
the plasma includes many particles with a small momentum perpendicular to 
the wall and so with a long wavelength perpendicular to the wall. When this 
wavelength is long compared with the particle mean free path, a "classical" 
treatment of the CP violation is not adequate. 

What follows is the description of a new method of calculation which apphes to 
all scenarios and all values of 1/ L. It reflects in the most direct way the interplay 
between the coherent phenomenon of CP violation and the incoherent nature of 
the plasma physics. This method can account for the generation of a CP violating 
observable from mass matrices with non-trivial space-dependence, as well as from 
particle interactions* In its simplest form, it easily reproduces the "thin wall" 
and "thick wall" calculations wih significant improvements over earlier estimates. 
In its more general form, it can incorporate effects which arise from the large 
diversity of scales present in a realistic plasma and can be the basis for Monte- 
Carlo simulations. For reasons of clarity, the method is best introduced with an 
example: the two higgs doublet model. The reader should bear in mind that it 
applies to other theories as well. 

7k- The latter mechanism dominates when the former one is not present. This is the case in the 
minimal standard model with CP violation originating from the quark Yukawa couplings.'"' 
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2. Construction of the CP violating sources 

Let us consider a set of particles with (not necessarily diagonal) mass matrix 
M{z) and moving, in the rest frame of the wall, with energy-momentum k. At 
their last scattering point Zo, these particles emerge from a thermal ensemble, prop- 
agate freely during a mean free time r ~ £, then rescatter and return to the local 
thermal ensemble in the plane Zg + tv, v being the velocity perpendicular to the 
wall, k±/E. During the time r, these particles evolve according to a set of Klein- 
Gordon, Dirac or Majorana equations coupled through the mass matrix M{z). It 
is in the course of this evolution that CP violation affects the distribution of these 
particles. Initially, at Zq, the contribution of these particles to any given charge 
cancel exactly the contribution of their antiparticles: (Q) = Tr[Q — Q] = 0; here, 
Q is the charge operator and the trace is taken over flavors as well as particle distri- 
butions. However, after evolving a time r across the CP violating space-dependent 
background, this cancellation no longer takes place for those charges which are 
explicitly violated by the mass matrix M{z). At the subsequent scattering point 
Zo + TV, these charges become (Q) — Tr[A^QA — A^QA] and assume a non-zero 
value, as A, the evolution operator over the distance tv, is distinct from its CP 
conjugate A. 

To be speciflc, let us deflne J±, the average current resulting from particles 
moving toward positive(negative) z between Zq and Zq+A, A = tv. The current J+ 
receives contributions from either particles originating from the thermal ensemble 
at point Zo, moving with a positive velocity and being transmitted at Zg + A, or 
from particles originating at Zg + A, moving with velocity —v and being reflected 
back towards Zq + A(Fig. la). A similar deflnition exists for J_(Fig. lb). J± are 
CP violating currents which are associated with each layer of thickness A moving 
along with the wall. Once boosted to the plasma frame, these currents provide CP 
violating sources, which fuel electroweak baryogenesis. 

The calculation of these currents is facihtated by the use of CPT symmetry and 
unitarity. CPT symmetry identifles the amplitude for a particle to be transmitted 
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from the left with the amphtude for its CP conjugate to be transmitted from 
the right, while unitarity relates transmission to reflection amplitudes. Instead of 
writing a cumbersome but general formula for these currents, let us work them out 
for a specific situation: a single fermion with a Dirac mass M{z) = m{z)e'^^'\ It 
could be a top quark having its mass generated from a two higgs-doublet lagrangian 
with an explicit CP violating term in the higgs potential, in which case, r is the 
mean free time for quark-gluon scatterings. As for the current, we choose the axial 
current ^ For this situation, the four- vectors J± take the form 




In this expression, w , = k_\_/E, is the velocity perpendicular to the wall at point 
Zo, V is the velocity a distance A away, = + (m^(zo) — m?'{zo + A))/E'^ and 
n{E,v) is the Fermi-Dirac distribution nj boosted to the rest frame of the wall, 
= [exp[7^£^(l — vvw)] + Q{zo,k,T) is the charge asymmetry which results 

— * 

from the propagation of particles of momentum k in the interval [zq, Zq + A]. In 
our specific example, Q is the chiral charge and is given by 

Qa{zo, k, t) = |T^|2 - \Tr\^ - |T;^|2 + |r^|2, (2) 

where Tl is the amplitude for a left-handed spinor to propagate over the distance 
A, Tr = Tl{M -Mt) and Tj^ = Tl{M M*). Finally, the brackets (. . 
in Eq. (1), average the location of point Zg within a given layer of thickness A as 
scattering occurs anywhere within a layer. 

t We choose this current because in the two Higgs model a combination of axial top num- 
ber and Higgs number diffuses efficiently into the symmetric phase and is approximately 
conserved by scattering in the symmetric phase.'*'^^' 
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The standard "thin wall" and "thick wall" situations are obtained in taking 
t/L to cxo and 0, respectively. In the "thin wall" limit, r/L — > oo, the amplitudes 
T become the usual transmission coefficients and our expressions (1) for J±, match 
trivially earlier calculations of scattering of particles off a sharp interface. In the 
"thick wall" hmit, r/L — > 0, the currents J± yield, after a boost to the thermal 
frame, a locally defined space-time dependent source density S{x, t) which gener- 
alizes, and gives a precise meaning to, the local CP violating operators already 
considered in the literature. In our example, the source per unit volume per unit 
of time, located at a point x fixed in the plasma, at any given time is, to first 
order in vw and r/L, 



1 Zo+Tv/2 OO ^ 



Zo—Tv/2 'ym{z) 



Zo=^w{x—Vn!t) 

(3) 

In order to obtain an explicit form for the source Sa{x, t), we need to compute the 
CP violating charge Qa and perform the integration in Eq. (3). 

3. Computation of the charge asymmetry Q 

In general,Q(^o, k, r) is a charge asymmetry produced by particles moving with 
momentum k between the planes Zq and Zq + tv. Its calculation may require a 
different technique depending on the relative values of the time scales involved 
and on the choice of the charge. The physics of the generation of a CP violating 
observable, is the physics of quantum interference. It is most easily dealt with by 
treating the mass M{z) as a small perturbation {M{z) < T). Using techniques 
developed in Ref. 4, one finds, for the transmitted amplitude, 

^o+A Zl 

Tl{zo,t) = e'^''^ [1 - J dzij dz2 e^'^k^^^^-^^) M{z2)M\zi) + 0{M/k^_f]. 

(4) 

This expression has a straightforward diagrammatic formulation presented in Fig. 
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2a. A similar expansion can be written for the reflection amplitude Ri (Fig. 2b). 
The various terms in the sum correspond to various paths with different CP odd 
and CP even phases. Only interference between these paths contributes to a CP 
violating physical observable such as Qa{zo, k,T); one finds 

Zo + A Zl 

QAizo,k,T) = 8 J dzi J dz2 siYi2k^{zi-Z2) \m[M{z2)M\zi)] + 0{M/k^f. 

Zo Z 

(5) 

This expression is valid for any wall size and shape, and generalizes easily to the 
case where many flavors mix. Given a wall profile, the integrals can be evaluated 
and Eq. (5) can be inserted in formulas (1) to provide an explicit form for the 
currents J±. 

It is simplest to work out the case of a very thick wall L ^ t. Using the 
derivative expansion M{zi) — M{zo) + (zq — Zi)dzM{zo), one finds 

Qa{zo, k,T)^- Af{k^A)/kl Im[Mt 9,M],„ 

= - Af{k^A)/kl m2 d,9U^ with /(O = sine (sin^ - ^ cosC). 

(6) 

Inserting this latter expression into our formula (3) for the source density Sa yields 

SA{x,t) = -T^^,Vu,m''d,Jl^^^^^^_^^^^ X ^I{T,m,T)+o(^vl,{m/T)^{T/Lfy 

(7) 

X(r, m, T) is a form factor whose general form is 



rT(j/-(f)2i) 

If f dt 

I{r,m,T) I dy ^ ^-^^ J fit) 

m 

T 

with = + M|, . 

We have included thermal corrections. My, in the mass dependence of Eq. (8)to 
take into account modifications of particle dispersion relations from scattering. The 
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effects of scattering on particle propagation can be accounted for by substituting 
quasiparticles for particles, in which case, r is to be replaced with I/27, where 
7 is the width of the quasiparticle. Correspondingly, the dispersion relation is to 
be modified to incorporate self-energy thermal corrections. In the particular case 
of quarks scattering off gluons, the width 7 is ~ QsT/S, while the main thermal 
corrections amount to the shift ^ + M|,, with M|, = ggT^/B ~ 

The form factor X is plotted in Fig. 3. X vanishes as r ^ 0, it peaks at tT ~ 1 
and is well approximated by ~ 1/VtT in the range tT > 5. The interpretation of 
this behavior is straightforward. As explained earlier, constructive interference is 
maximal for particles whose transverse Compton wavelength kj^ is of the order of 
T, that explains the peak at tT ~ 1. As tT increases, fast oscillations along the 
distinct paths tend to cancel against each other and the resulting asymmetry drops; 
as a matter of fact, in the extreme limit tT ^ 1, the propagation is semi-classical 
and the asymmetry vanishes as it should.^ 

In the opposite limit, as tT 0, the asymmetry vanishes as the quantum 
coherence required is washed away by the rapid plasma interactions. Fig. 3b 
demonstrates the mild dependence of X on in. 

For the sake of comparison, we present an approximate form of Eq. (7), valid 
for m <C T, 

SA{x,t) -i.^=jn^Vn,T m^d.e for tT > 5 (9) 

SA{x,t) ^ --^ ^^VyjT m^dz9 for rT - 1 - 2. (10) 



* A systematic method which accounts for the thermal structure of a quasiparticle in the 
interference mechanism is given in Ref. 4. 

f However, a semi-classical treatment alone might miss the important contributions of long- 
wavelength particles moving at large angles in respect to the wall motion; without their 
contribution, the asymmetry would fall off as fast as 1 /r. 
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4. Application to two-higgs models 



The formula we derived for the chiral source Sa{x, t) can be directly applied to 
the top quark propagating in the thick wall of a bubble produced at the electroweak 
phase transition in two-higgs models. Here the mean free path is dominated by 
gluon scattering r ~ 3/(2ggT) ~ (1 — 2)/T and is typically smaller than the 
estimated thickness of the wall: t/L 0.01 — O.l.'^'^'The wall velocity is jw^w < 
l!''"'" Finally, The mass of the top quark mtiz) is Yt(f) <YtTc^T while My ~ T/2, 



hence, m = y m| + Mj, < I.IT and the assumptions under which we derived Eq. 
(10) are approximately fulfilled. We find 



where the number of colors Nc — S. 

Recent work on the source terms for axial top number in the two Higgs model in 
the thick wall case, have treated the fermion interaction terms with the background 
Higgs field as a CP violating contribution to a classical Hamiltonian in computing 
the CP-violating perturbation to particle distributions. Ref. 8 pointed out that 
these interactions split the energy levels of particles and their CP conjugates in a 
way reminiscent of a chemical potential. These classical treatments obscure the 
origin of the CP violating effect as resulting from quantum interference. However, 
these methods, if implemented properly, should provide reasonable approximations 
to our formulae for those particles whose wavelength l/k± is short compared with 
VT. As an illustration, Ref. 13 found a source term SA{x,t) = 1/3 T VwUi^dzO. 
These authors did not account for the quark-gluon interactions in the rate for 
incoherent axial top number violation ( a factor 1/r) and the three dimensional 
phase space (a factor of ~ 9({3)/'k'^), factors which arc numerically unimportant 
but which are needed for theoretical self- consistency. Even after including these 
effects, our formulae do not agree for large r because our integral (3) is dominated 
by particles with long wavelengths in the direction perpendicular to the wall, for 




SA{x,t) 




(11) 
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which a classical approximation is not adequate. Numerically, for tT ~ 2, our 
answer approximately agrees in magnitude. 

5. Conclusion 

In conclusion, we have introduced a new method to compute CP violating 
sources resulting from particle interaction with an expanding bubble during a first- 
order electroweak phase transition. 

1. The method refers to explicit physical processes in the plasma. In particular, 
it does not make use of thermally averaged operators, or effective chemical 
potentials whose connection to the microphysics is indirect, as they do not 
vanish as the mass m and/or the mean free time r vanish, and whose appli- 
cability is restricted to the range L ^ t ^ 1/T. 

2. The method makes explicit the quantum physics of CP violation and its 
suppression resulting from thermal effects. 

3. A major advantage of our formulation is that it easily applies to charges gen- 
erated by flavor mixing through arbitrary large mass matrices. In particular, 

[17] 

it can be applied to cases, such as the supersymmetric standard model for 
which there is no known semi-classical approximation. 

4. It is vahd for all wall shapes and sizes as well as for arbitrary particle species 
and interactions. 

5. Finally, it incorporates CP violation which originates from particle interac- 
tions as well as from non-trivial space-time mass dependences. In particular, 
it generalizes and agrees with the decoherence arguments invoked to rule out 
electroweak baryogenesis from CP violation in the quark Yukawa couplings, 
as given in Ref. 4. 
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#DE-FG06-91-ER40614. The work of A. N. was supported in part by a fellowship 
from the Sloan Foundation. 
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FIGURE CAPTIONS 

1) (a) Contributions to J+. (b) Contributions to J_. 

2) (a) Amplitudes contributing to T^. (b) Amplitudes contributing to R^. 

3) (a) The form factor X is plotted versus rT, for the case m <^ 0, My = 
T/2 (quark-gluon interactions). The dotted curve results from numerical 
integration of Eq. (8); for values of tT > a few, it is well-approximated by 
its asymptotic form 1/2\/tT (solid line), (b) The dependence of X on the 
mass m is mild in the range m < T. 
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Fig. 2 



I 



0.35 
0.3 

0.25 
0.2 

0.15 
0.1 

0.05 



(a) 



l5^-=0.5 
T 



10 



15 



20 T T 



0.25 



0.2 ■ 



0.15 



0.1 



0.05 




0.2 



0.4 



0.6 



. 8 



T T=10 



m 
T 



16 



Fig. 3 



